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Abstract 

Products of random transfer matrices are applied to low dimensional disor- 
dered systems to evaluate numerically extensive quantities such as entropy 
and overlap probability distribution. The main advantage is the possibility 
to avoid numerical differentiation. The method works for arbitrary disorder 

distributions at any temperature. 
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Products of random transfer matrices are a powerful tool for the study of low dimensional 
systems Up to now they have been mainly used to evaluate free energies and correla- 
tion functions. From a formal point of view, the other thermodynamical quantities can be 
obtained by differentiating the free energy with respect to suitable parameters. However, in 
most of the cases, one has only numerical solutions and, expecially at low temperatures, a 
very high precision is required for a precise differentiation. In practice, it is very difficult 
to estimate quantities such as the entropy or the specific heat even in a one dimensional 
disordered system. 

The purpose of our letter is to introduce a new numerical method to evaluate ther- 
modynamical quantities directly from the products of random transfer matrices avoiding 
numerical differentiation. The method is motivated by the papers of Masui et al. p|J^ on 
the study of the number of metastable states in one and two dimensional Ising disordered 
systems. As specific example we shall discuss the overlap probability distribution P{q) for 
zero and finite temperature for a one dimensional Ising chain. The method can be extended 
to higher dimensional systems with short range interactions. 

An Ising chain consists of N Ising spins with nearest neighbour interactions in an external 
field. The interactions and/or the local field are random quenched variables. At zero tem- 
peratures the only states which are important are the metastable states in which the spin cTj 
is aligned with the total magnetic field hi acting on it. If we denote by the spin configura- 
tion in the metastable state a, the overlap between pairs of states is Qap = J2i 
where a and (3 denote a pair of metastable states, and the sum is extended to all the 
spins. In general g^/j depends on the realization of disorder and on the pair of states. The 
number of metastable states with a given overlap can be obtained from 

1,N 1,2 

AA^(^,/5) =EEe^^''" n n ^(^>r) (1) 

where and are the hamiltonian of two identical replicas of the system. The step 6- 

functions ensures that only the metastable states are counted. The Gibbs weight exp(— 

has been introduced to select the energy of the metastable states. For any given value of (3 
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the sum Eq. (|l|) will be dominated by metastable states of defined energy. In the same way 
any given value of uj selects metastable states of defined overlap. Due to the 9 function in 
Eq. (|I]), the parameter [3 is not necessary related to the temperature, but it can be seen as 
a Lagrange multiplier. 

N'n{'^i the number of metastable states of a chain of length N for a given value of j3 
and uo, can be written as a sum over the energies and overlaps 

Un{u:, /3)=Y: ^fNiq, E) e-^^^^^-\ (2) 

For large N the sum is dominated by the saddle point energy E' and overlap g', so that we 
can write 

Ar^(c.,/3) ~ Ar^(g',E')e-'^^'+^-'^'. (3) 

Equation (|l]) can be written as a product of random matrices, so that the Oseledec theorem 
ensures that in the thermodynamic limit the quantity in Eq. is the same for almost all 
realizations of disorder Consequently one has the Legendre transform 

\^M^{q, e) = ^ In Ar^(^, /5) + 2 /? e - a; g (4) 

where the bar denotes the average over the realizations of disorder. Here e = eiio^fi') is the 
energy per spin and q = q{uj,(3) the overlap selected by the chosen value of uj and [3. In 
particular the average overlap is obtained as 

q{uj = Q,l3) = q = Jqmdq (5) 

and is the value of q where limN^oo{^/N)\nNN{q,e) reaches its maximum value. 

The key point is that not only {l/N)\iaAfN{^^, P) but also q{uj,P) and e{u,l3) can be 
computed directly from a product of transfer matrices. 

We shall denote by Si = {crj, af) the pair of spins at the same site in two replicas of the 
system. Equation (|^) can then be written as a chain sum over Sj, i.e. as a product of the 
matrix {Sij^i\Ni\Si) . From Eq. (|l]) it is easy to see that iVj obeys the recursion relation 
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Si-i 

X Ee{hyr)e-'^" (6) 



with an initial matrix A'^^ with all elements equal and with their sum equal to unity. Here 
Ef is the energy of the spin erf. The matrix A^j'_i is obtained by rescaling A'j.i by the factor 
nj_i given by the sum of all its elements. This ensures that the sum of the elements of A^j'_i 
is one. 

The matrix A'^j gives the number of metastable states of a chain of i spins for all the 
possible configurations of Si and Sj+i. Thus the logarithm of the number of metastable 
states per spin, averaged over the disorder, is 

1 1 N 

-lnAA^(a;,/?) = Jim -^Inn, (7) 

1=1 

The value of q and e is obtained from the derivative of Af]y{u, (3) with respect to u and /3, 
respectively. These can be expressed in terms of the matrices derivative of iVj. For example, 
defining the matrix Oi as {Si^i\Oi\Si) = (d/duj) (S',;_|_i|iV'j|S'j) the derivative of Eq. (|]) yields 
recursion relation 



(5,+i|C),|5,) = ^ [(5,|6:„i|5,_i) 

+aV2(5,|iV;_i|5,_i) 



xe!^'^i-l n 0{hy^)e-^''^ (8) 

o=l,2 



with the initial condition 

{S2\d,\Si) = crlal{S2mSi). (9) 
The matrix 0^_^ is related to Oj_i as 

{S,+M\Si) = [(5,+i|6i|5,) - o,(5,+i|iV;|5,)] /m (10) 
where Oj is the sum of the elements of Oj. 
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From the renormalization Eq. ([lOD follows that the value of q for the chosen value of uj 
and /5 is given by the weighted average of Oj as 

1 ^ Oi 

q = lim — —. (11) 

1=1 ' 

A similar equation yields e. 

The method can be generalized to finite temperatures. In this case we have to count all 
the states, and not just the metastable ones. Thus We should thus eliminate the 6'-functions 
in the above equations. Moreover, in this case f3~^ is the temperature of the system. 

In the absence of the ^-functions, which relates the states in z — 1, i and z + 1, we can use 
a simpler algorithm 0. The equations are formally the same, but with matrices replaced 
by vectors. Therefore, in analogy with the zero temperature case, we introduce the vector 
Ni{Si^i) which now obeys the recursion relation 

iV.(^.+i) = E e-^^^^^^-^^'^? N'US,) (12) 

where the initial vector N'q{Si) has all the elements equal to one. At each iteration the 
vector is normalized so that the sum of the elements of N[_-^^{Si) is one. To evaluate q 
and e we introduce the vectors of the derivatives of A/"j(S'j+i) with respect to u and /3, 



respectively. For example, taking the derivative with respect to oo of Eq. ([121) leads to the 
vector Oi{Si+i) = {d / du) N i{S i+i) obeying the recursion relation 

^e-''>" ne-^^" (13) 

a 

with initial condition 0'q{S\) = 0. The vector Oi and O- are related by 

0:(5,+i) = [0,(5,+i) - o, N[{S,+,)] Im (14) 

where as above rii and Oi are the sum of the elements of Ni{Si+i) and Oi(S'j+i), respectively. 
The average energy e is evaluated in a similar way. 
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The specific heat and susceptibihty can also be computed with this method by introduc- 
ing the matrix or the vector of second derivative. 

We have analyzed only Ising chains, but our arguments apply to any formally one di- 
mensional system and can be easily generalized to two and three dimensional systems by 
considering strips and bars, respectively. The numerical calculation becomes, however, very 
heavy since the dimension of the transfer matrices and vectors grows exponentially with the 
transversal dimension of the system. 

Let us apply this method to the Ising chain in a random field described by the hamiltonian 

N N 

H = - ^Jaiffi+i-^hiai (15) 

i=l i=l 

where hi is a random local field. Without loosing in generality we can take J = 1. It is 
known that if hi are distributed with probability distribution 

V{hi)=p6{K-h) + {l-p)8{K + h) (16) 

with < p < 1 then the zero temperature entropy is different from zero. This is peculiar of 
distribution (16). In fact, if we replace the delta functions by two sharp peaks of width A/i, 
then S{T) oc T for T ^ A/i This drastic change can be clearly seen in the behaviour of 
q = q{uj = 0, f3 = l/T) as a function of T, as shown in Fig. |l|. In Fig. ^ we report the spin 
glass susceptibility xsg, defined as the derivative of q{uj, (3) with respect to uj ioi u = 0, as a 
function of temperature for two different values of A/i. The value of xsg has been obtained 
by using the vector of the second derivative of Ni with respect to uo. Finally Fig. |^ shows 
the typical form of limiv_»oo(l/A^)ln A^Ar(g, e). The maximum is attained for q = q. In the 
figures J = 1 and h = 2. The convergence of the method is very fast. Reasonable good 
numbers are obtained already with O(IO^) iterations. For T = 0, the results are in very 
good agreement with the theoretical prediction of Ref. 0. 

The parabolic shape of lim7v->oo(l/A^)ln A^Ar(g, e) means that the overlap probability dis- 
tribution P{q) is selfaveraging and equal to a ^-function on the average value q. This is a 
general result for one dimensional systems with short range interactions 0. 
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We conclude by noting that since for T — the probabihty distribution P{q) remains 
a 5-function the hmit T ^ is not singular and the results of the zero temperature 
method perfectly reproduce those of the T — limit of the finite temperature method. 

In conclusion we have described a method to evaluate all thermodynamical quantities 
in terms of products of random transfer matrices avoiding numerical differentiation. The 
method can be used, in principle, also for two and three dimensional systems. The method 
is not restricted to spin systems, but it can be used in any other problems which involves 
products of transfer matrices, e.g. directed polymers. 
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FIGURES 

FIG. 1. The average overlap g as a function of temperature T for = (a) and = 0.3 
(b). In both cases J = 1 and h = 2. 

FIG. 2. The spin glass susceptibility xsG as a function of temperature T for Ah = (a) and 
Ah = 0.3 (b). In both cases J = 1 and h = 2. 

FIG. 3. \imN^oo{l/N)\nNN{q, e) as a function of q for A^ = 0, T = and e = emin, the ground 
state energy. The maximum is obtained for q = q. In figure J = 1 and h = 2. 
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